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. Abstract. The specified constant 4- vector field reproducing the spherically symmetric 

I stationary metric of cold dark matter halo in the region of flat rotation curves results in a 

^ ■ constant angle of light deflection at small impact distances. The effective deflecting mass 

OO , is factor 7r/2 greater than the dark matter mass. The perturbation of deflection picture 

due to the halo edge is evaluated. 
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1. Introduction 

For spherically symmetric stationary dark halos in spiral galaxies one could explore the 
■ weak field approximation of pressureless matter implying the metric 

.F^.: d.^ = {1 + 2 </>GR(r)}dt2 2GMuu{r) ~ "'f^^' + ''''' 

r 

where MdmI^) denotes the dark matter mass inside the sphere of radius r, since to linear 
order in small MDM(yr)/r and 0gr(^) the Einstein equations give 



SnGpir) = G\ = R\-\ g\ R = ^-^^^M (2) 
exactly yielding 



r 

J dVp(r) = MDM(r), (3) 



while 0gr(^) represents the gravitational potential of dark mass 
... ? , GMDM(r) 

(i)GR{r) = dr — (4) 
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producing the force 

F = V</)GR(r) = (5) 

as well as a negligible pressure 

87TGpn{r) = -Gl^0, n = {r,e,y^}. (6) 
For fiat rotational curves with the characteristic velocity Vq we get 

vi GM^uir) 



(7) 



and, hence, 



G 

which argues for the metric 



MuMir) = ^vlr, (8) 



ds' = f (r) dt' - dr' - r' [dO' + sin' 6 dip'] , (9) 



with 

2v' 

h(r) = 1 - 2vl f'(r) = (10) 



We consider the metric is close to the Minkowskian one up to small corrections in Vq 0, 
so that 

f(r) = 1 + 2t;olnr. (11) 

Surprisingly, metric was derived for the dark matter modelled by 4-vector field Au 
with constant components [1] 

A = ^^, A=^^-. (12) 

In [2] we supposed a treating the vector field in terms of ghost condensate [3-10]: dt(f) = 
(f) = Ao, and global spherically symmetric monopole il^{r) = Vx('") = -^if) ^ n = r/r, 
combined into the 4-vector A^- The phenomenon of ghost condensate is an analogue of 
Higgs mechanism for the condensation of scalar field with negative square of mass near 
zero values of field. This mechanism of ghost condensation usefully works in the case 
of expanding homogeneous universe, that has a specific reference frame with a selected 
temporal component in comparison with spatial ones. Then, the scalar field with negative 
sign of kinetic term near zero values of kinetic energy in its action composed of invariant 
derivative terms X = c?^0(5'^0, condenses at a nonzero value of <p, which represents a 
temporal component of 4-vector field. The ghost condensation removes the problem with 
negative kinetic energy, since near the condensate the fiuctuations have correct sign of 
kinetic energy, of course. The studies [1] show, that field 4> should be constant due to the 
current evolution of universe, with a high accuracy. Next, one could combine the ghost 
condensate with the global monopole, that produces a correct static spherically symmetric 
profile of mass distribution p ~ 1/r^ in the region of fiat rotation curves due to the dark 
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matter. Therefore, one should suppose the presence of global monopole in the center of 
spiral galaxy, i.e. inside a supermassive black hole. At large distances, the monopole 
corresponds to the spatial component of vector field. Then, one considers the background 
field (f) = fi^t — fi^r at large distances (perturbations become important at small distances, 
corresponding to inner edge of halo), where /x* is the scale about the Planck mass, while 
the ratio of determines the rotation velocity Vq (see details in [2]). So, the ghost 

condensate in the presence of global monopole with account of higher derivative terms in 
the action can be represented by the constant 4- vector field relevant to the situation of 
both universe expansion and dark halo. 

In the present paper we calculate the light bending by the metric modelling the dark 
halo. In Section 2 we in detail calculate the defiection angle in the limit of 'infinite' 
halo implying a negligible value of impact distance with respect to the halo size and 
compare the result with the deflection by a global monopole [11]. Then we determine 
the effective deflecting mass and flnd its difference by a factor of 7r/2 from the dark matter 
mass responsible for the flat rotation curves in contrast to naive expectations of their 
coincidence. Further, we describe the geometry of deflection revealing properties consistent 
with astronomical observations. In Section 3 we calculate the deflection in the case of flnite 
halo, that demonstrates a weak dependence on the impact distance slowly perturbing the 
geometry discussed in Section 2. The results are summarized in Conclusion. 

2. Limit of infinite halo 

By the Hamilton-Jacobi equation 



incorporating two integrals of motion in the spherically symmetric static gravitational fleld 
at flxed polar angle = const., so that £ is the conserved energy of particle and Wl is its 
rotational momentum, which we put positive for deflniteness, we deduce 



g>'''d,Sd,S = 
for an action S of massless particle written in the form 

s ^ -st + me + s(r), 



(13) 



(14) 




(15) 



which results in 



r{t) 




(16) 



where V"^ 



is an analogue of centripetal potential. 
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Thus, the quantity S{r) represents the action for the radial motion of particle in the 
gravitational field with account of energy due to the angular rotation. The trajectory 
is implicitly determined by the equations^ 

r{t) 



dS 
dS 

dS 

dm 



const 



-t+ dr 



£ 



ro 
r{t) 



const = 9 



dr 



1 f m 



fh H _ ^2(^) 



Taking the derivative of p7|) and (fTBj) with respect to the time§, we get 

S 



1 
d 



r 



r 



fh y/£^-V^{r) 
r i m 

~2 ~ 



and, hence, 



^20 



(17) 



(19) 
(20) 

(21) 



relating the energy and rotational momentum. Let us denote the ratio of two integrals by 

m 

with the dimension of length, which we call as 'an impact distance', that becomes clear 
from Fig. shown for the fiat Minkowskian space-time. 
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Figure 1. Geometrical relations between the impact distance a and circular velocity r6. 



At the trajectory given by (fT^ and (pUj) we get 

d^ a If 1 
dr 



r^ \ h 



1 - 



(22) 



I Here we take the positive sign of square root for brevity of record. 
§ As usual dtf{t) - /. 
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where the negative sign corresponds to the branch of trajectory for the particle arriving 
from infinity, while the positive sign does to the branch for the particle departing to infinity. 
Denote the solution of 

= a^f (r) 

by '"min- At this value the radial velocity r becomes equal to zero (the denominator in 
Eq.((22I) vanishes), that corresponds to returning point. So, the light ray approaches the 
gravitating object from the infinity, while it reaches the distance rmin, after that the ray 
moves off the object. Therefore, the quantity rmin denotes the minimal distance between 
the light trajectory and the gravitational lens. Then, the overall change of 6 for the particle 
is given by 

a f 



Mg = 2 / dr — W , (23) 

J r2 V h / 

if we put the size of halo much greater than the impact distance: a^f/r^ — > 0, that 
corresponds to the limit of infinitely large halo. In contrast to [12], we expect that this 
limit is relevant to actual astronomical observations. 

Let us integrate out by expanding in small parameter of 0. So, with the 

accuracy up to the first order we get 

f = f (rmin) + 2vl In ^ + 0«), 

^min 

h = 1 - 2vl 
that allows us to write down 

Mg = AOo + v^oiAOi + + A^s) + C(wo), (24) 
where the leading term 

oo 

Aeo = 2 [ dr^v/f(r^i,) w ^ ^^^^ 



after the substitution 
a 



z=- Vf(r^in) (26) 



is reduced to 

1 



dz 

AOq = 2 , = 2 arcsin z 



(27) 











which corresponds to the straightforward motion with no deflection. The contribution A^i 
originates from the factor of 1 / ^/h ^ 1 + Vq in the integrand of (j^^ , so that 

A^i = A^o = TT. (28) 
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The expansion of 



Vf = V firkin) -v^lnz + Oiv^) 

results in 

1 



where we have neglected the deviation of f (rmin) from unit 

i{rnnn) = l + 0{vl). 

Analogously, the expansion of 

1 1 ^ z'^lnz 4x 

^0 ... + O{v^o) 



produces 



A^3 = -2 /d.^^il=. (30) 



7(1-^ 



2^3 



Summing up and (jHUI) gives 



A», + A03 = -2/d^-^=i=. (31) 



The integral is easily calculated by parts 



, In^; 2z f dz 

dz — , = , mz + 2 



2; 



2 ' 



so that 

A92 + A93 = 71. (32) 

Thus, the deflection of light is determined by the angle 

= AOg - tt = vliAOi + A^2 + A^g) = 2nvl. (33) 

Note, this limit of deflection angle is independent of impact distance. This result is similar 
to the calculation of light deflection by the global monopole metric in [11], where the radial 
component h can be taken in the same form, while the temporal component of metric 
is constant, f 1-^ const., which corresponds to the term of A^i in calculations above. 
Therefore, we emphasize that the deflection angle in the metric of global monopole is twice 
less than the calculated value of AO = 27[Vq. 

The bending angle derived is in agreement with results of [13], where the deflection due 
to the dark matter with a significant pressure was calculated numerically at various equation 
of state parameters m„ = Pn/ P, and the limit of pressureless matter was mentioned, too. 
The same result could be extracted from [12] by calculating the limit of f q — at infinitely 
small impact distances from analj^ical expressions in [12] also discussed in Section 3. 
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2.1. Deflecting dark matter 

The main feature of deflection angle under study is its independence of impact distance. 
Such the deflection can be ascribed to an effective dark matter mass. Indeed, let us compare 
the effect with the deflection by a nonrelativistic spherical mass distribution described by 
the Schwarzschild metric: 

4MG' 

A9bh = ■ (34) 



Putting 



T ■ 



we derive the expression for the effective dark matter mass responsible for the deflection 

MdmW^^^o'^, (36) 

which linearly depends on the distance. Compare it with the dark matter mass responsible 
for the centripetal acceleration in the region of flat rotational curves, 

— = 2 ^ AfoM r = TT^'o^- 37) 

Therefore, in the metric under consideration the effective deflecting mass of dark matter is 
factor tt/2 greater than the dark matter mass producing flat rotational curves in Newtonian 
limit. 

The same conclusion is drawn from virial estimates of dark matter mass in galaxy 
clusters, where we get the relation 



2, G'(Mdm) 



2(r) ' 

since the average kinetic energy of N galaxies with average masses (Mgaiaxy) = {M-qm) /N 
at > 1 

is related with the average potential energy of galaxies posed at average distance (r) 

^ > 2 (r) 2 (r) 

by the virial theorem 

2(7kin) + (V) = 0. 

Both estimates produce the similar linear dependence of dark matter mass versus the 
distance, as we retrieve from the gravitational lensing. 

Remarkably, astronomical observations of galaxy clusters give 

(Mdm) ~ 200(Mb), (38) 
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where (Mb) is the mass of luminous baryonic matter in the clusters. Note, that the mass 
of luminous matter is related with the velocity of rotation in spiral galaxies at distances of 
flat rotational curves by the empirical Tully-Fisher scaling law 

(Mb) ~ ^;o^ (39) 

which is natural, if one suggests the existence of universal critical Milgrom acceleration [14] 

ao^Hi.£(^, (40) 

where ro is a characteristic distance for a galaxy, where the dominance of baryonic matter in 
the gravitation comes to competition with the dark matter contribution becoming dominant 
at larger distances. The value of tq is close to the visible size of galaxy, i.e. it is the 
characteristic size of baryon distribution in the galaxy. So, let us compare the masses of 
luminous and dark matters: 

(Mb) ~ ^ vl ro, (Mdm) - ^ vl r„ (41) 

where r^ is the size of dark halo. Therefore, we conclude that 

r-^~ 200ro, (42) 

i.e. the dark halo spreads far away from the visible size of luminous object. This point 
supports the relevance of considering the limit of infinite halo in the problem of gravitational 
lensing at realistic impact distances about the visible size of lens. 

Thus, the vector field model of dark matter qualitatively reproduces both dark matter 
effects for astronomical objects: the fiat rotational curves and light bending. The numerical 
effect on about 50%-difference between the effective deflecting mass and dark matter 
gravitational mass in rotational curves could probably be used in future for falsifying the 
model, though such a combination of astronomical data seems to be quite exclusive. 

Another aspect of comparing the rotation curve picture in a galaxy with the 
gravitational lensing by the same galaxy has been recently considered in [15], where the 
potentials responsible for the rotation and lensing have been introduced. The difference 
between the potentials would signalize on a non-zeroth pressure of dark matter. In the 
case under study these potentials coincide, so the lensing of dark mass distribution could 
be evaluated by standard techniques for the cold pressureless matter. In contrast to the 
potentials, we operate by the effective deflecting mass allowing us to compare it with the 
mass one could expect from the rotation velocity. This could be important for a treatment 
of some misalignment between the lensing mass and dark matter mass in astronomical 
measurements of X-band lenses [16]. 

2.2. Geometry of deflection 

Let us consider the geometric aspects of gravitational lensing. In Fig. |2l point S denotes 
a star, which light reaches an observer posed in point O by path SBO curved by dark 
matter object somewhere placed at interval BD. Since all of angles are extremely small, 
while distances between the observer, dark matter object and star are astronomically large, 
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we can perform calculations in the first order versus angles. Then, the distance between 
the deflection point and position of dark matter object projected to the line between the 
observer and star r = BD (see Fig. ^ equals 

r = la] r = d(3. (43) 

Points D and D' coincide because \DD'\ I = 0{a'^) is negligible. The angle of 
deflection Ad is independent of r and 

M = a + p. (44) 

Therefore, 

/J = M^. (45) 
The same deflection takes place by the path in the upper half of plane shown in Fig. |21 




Figure 2. The plane geometry of light bending. 




Figure 3. The gravitational lensing of linear image by the dark mater object {left) in 
comparison with the initial pattern without the light bending [right). 

Thus, the angle between two images of star is given by 

A(3 = 2(3 = 2 Ae (46) 
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Figure 4. The gravitational lensing of centered circular image by the dark mater object 
(left) in comparison with the initial pattern without the light bending (right). 

This geometry can be easily applied to draw samples of gravitational lensing in the model 
with constant A6. 

So, Fig. El represents the gravitational magnifying of linear image at c? = 1, / = 10 and 
small angles. The picture illustrates the characteristic distortion caused by such lenses. For 
simplicity of picturing, in Figs. 3-6 we put the parameter Vq = l/(407r) ^ 0.008 ~ (0.09)^ 
in order to enlarge the visual effect. Nevertheless, in practice, an angular magnification of 
real astronomical facilities is sufficiently large for the observation of gravitational lensing 
by the dark matter (see numerical estimates in the end of section). 

One can easily recognize basic features of deflection and distortion: 

(i) The dark matter object is posed in the center of inversion symmetry. 

(ii) Hence, the angle distances between the points inside two images coincide. (If one 
displaces the star image by angle 7, then the second image is displayed to the same 
angle increment 7.) 

(iii) Therefore, small radial widths of images are coincide, while the radial sizes are 
proportional. The coefficient of proportionality k is equal to the ratio of angle fractions 
given by division of A/3 by the center of symmetry, i.e. the dark matter object. This 
point is shown in Fig. |3 

(iv) Then, the ratio of surface areas is equal to k as illustrated in Fig. El for the distortion 
of circle object. 




Figure 5. The gravitational lensing of circular image by the dark mater object (left) in 
comparison with the initial pattern without the light bending (right). 
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Thus, for the model of gravitational lensing, we arrive to the following conclusion: 
The relative values of surface brightness for two images of lensed small object is given 
the inverse ratio of angle fractions, k, or visible impact parameters. 

So, the ordinary picture of gravitational lensing due to the dark matter modelled by 
the vector field is shown in Fig. IHl We have taken two galaxies as elliptic figures at different 
distances (ci = 1, Zi = 9, /2 = 3). Their images are narrow arcs surrounding the dark matter 
object, that is generically similar to the real astronomical picture inserted as illustration. 




• ■ A • s 

Figure 6. The gravitational lensing of sample galaxies at different distances by the dark 
mater object (top-left) in comparison with the initial pattern without the light bending 
(top-right) and with the real picture of Abell 2218 cluster lensing (bottom) [from NASA, 
1998]. 

2.3. Numerical values 

The characteristic angle between two images for an object due to the gravitational lensing 
is about 

Af3 ^ 2A9 = 4nvl ^ 2.6 ■ 10*^ arcsec. (47) 
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A typical velocity of rotation is in the range of 100 — 400 km/sec, i.e. Vq ~ 10 ^. Then, 

A/3 ~ 2 arcsec. (48) 

Such the value is consistent with astronomical observations of gravitational lensing, where 
the characteristic angle sizes between the dark matter object and image of lensed luminous 
object is about 1 arcsec. 

3. Edge of halo 

If the impact distance is comparable with the size of dark halo, then one has to take into 
account the halo edge at r = (see Fig. [7j). So, the formulae exploited above give the 
angle increment in the halo, 

(Z 111 z \ 
27r — 4 arcsin + 2 * * I , (49) 

with = a/r^. 



AO 



Figure 7. The propagation of light through the dark halo of finite size r^,. 

Let us suppose that beyond the halo we get the Schwarzschild metric with the 
gravitational radius r^: f,h i-^ 1 — y. Then, the increment of angle beyond the halo 
is determined by the integral 

oo 

A^BH(r.) = 2 / dr ^ , / (50) 

1 (l 
V V r 

which can be easily evaluated after the substitution 

1 
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and expansion to the first order in Vg/a. So, we get 



a 1 + ^ 



2 arcsin z-^, + 2 — { 1 — \/l — zl 
a 



13 



(51) 



The value of gravitational radius could be fixed by the dark matter mass at r^, 

rg = 2GMDM(r.) = 2vlr,. 
Then, the overall deflection angle 

A^, = A^G(r.)+A0BH(r.)-7r 

is given by 

9 I Zi,\n Zi, A c I 

Ae^ = vU2n-A arcsin z^ + 2 + — ( 1 - ^1 - 



(52) 



(53) 



zt V 

valid at ^ 2;^ ^ 1, while at 2:* > 1 we get the ordinary angle of bending A^bh = 2rg/a. 



Adjvl 



Adjvl 




0.25 0.5 0.75 




0.25 0.5 0.75 



Figure 8. The deflection angle light by the dark matter with the edge at r = r^, 
versus ~ ajvi, under matching with the Schwarzschild metric at Vg = 2vQri, (left) 



and 



3v^ri,/2 (right) in comparison with the deflection by the massive point of 



M = rg/{2G) [dotted line). 



Note, that the derivative of deflection angle has a singularity at 



dA^^ 



dz* Jl - zl V^o^vc 



2r„ 



1. 



(54) 



The singularity reveals as a fracture point of curve in Fig. |S| The fracture disappears 
at rg = 1^0?"^. However, such the redefinition of parameter r^, removing the fracture of 
bending curve (the derivative of A6'^ with respect to 2;^ becomes a continuous curve), results 
in a discontinuity of metric coefficient h at the edge of halo, that is absolutely unacceptable. 
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The origin of singularity in Eq. (|54j) is the specific dependence of halo metric-coefficient 
f on the distance r. So, it produces a slow deflection with negative sign in vicinity of halo 
edge {z^ = 1), i.e. there is a small repulsive component of force acting to the light ray in 
comparison with the case of point-like mass. It is clearly seen in derivation of Eqs. ()3Up 
and (jn^- Then, the attraction by the halo near the edge is small in comparison with the 
case of point-like mass, hence, the substitution of trajectory passing the halo instead of the 
same path in the field of point-like mass results in an essentially smaller deflection, and we 
even get the fracture of curve for the bending angle, while the deflection increases when the 
path in the halo becomes longer. Thus, the fracture is the direct consequence of metric, 
specifled for the dark halo. Moreover, one could speak about a gravitational refraction of 
light on the border of two media with different metrics (in the case under study, we get the 
transition from the vacuum to the halo medium). 

The analysis shows that in quite a wide region of impact factors, the deflection by the 
dark halo with the edge is slowly perturbed due to the edge: the dependence on the impact 
distance is weak in comparison with the rapid change in the case of point-like mass. So, 
the picture of deflection considered in Section 2 qualitatively remains the same. 

Comparison with the paper by Nucamendi, Salgado, Sudarsky [12]. In [12] authors in part 
considered the same problem of light bending by the dark halo with the edge. The difference 
is twofold. First, the metric coefficient f in [12] has been taken in the form 

f = A j . (55) 

Second, the deflection by the Schwarzschild metric calculated by expanding in powers of 
inverse distance has been written in a form, which is literally different, but in fact it should 
coincides with the formula (jHH) with the accuracy up to O^Vg/a"^). These facts allow us to 
make a direct comparison in the limit of — 0. So, we flnd that our result is consistent 
with the corresponding limit of expression in [12] with the accuracy mentioned above. In 
part, the limit of inflnite halo — > in [12], reproduces the value of = 2'kvq, which 
is actually in the evident contradiction with the main result of [12] shown graphically as 
the function of z^, (see Fig. 1 in [12]). Indeed, instead of 27r one could observe a number, 
which is factor 4/3 greater. The same factor remains, when we compare the complete curve 
in [12]: the deflection angle of dark halo with the edge is multiplied by 4/3. Literally, the 
curve of [12] could be reproduced from our right curve in Fig. |H1 multiplied by 4/3, that 
gives the correct value for the bending by the Schwarzschild metric, but it misleads on the 
deflection by the halo. The factor 4/3 removes the fracture point from the curve, of course, 
as it is clearly seen in Fig. 1 in [12]. 

Thus, our analysis on the deflection angle by the dark halo with edge is consistent with 
the formal limit of expressions in [12], but we disagree with the numerical in-flgure-shown 
result of [12]. 
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An inner edge. We have considered the external edge of halo at large distances. However, 
the luminous matter essentially contributes to the deflection at short impact distances, 
where the gravitational field is dominantly determined by distributions of stars and gas. 
So, there is the inner edge of halo. This problem has been comprehensively investigated in 
terms of Modified Newtonian Dynamics (MOND) [14]. Recently, Bekenstein formulated a 
consistent relativistic theory of gravitation including gravitational scalar and vector fields 
in addition to the ordinary tensor of metric [17], so that it incorporates general MOND 
effects originated from the introduction of critical acceleration oq. The authors of [18] have 
found the same value of deflection for the case of inflnite halo at impact parameters greater 
than the inner edge, that corresponds to the isothermal halo of circular velocity vq. The 
comparison with observational data has been performed, so that the conclusion drawn is 
the following: MOND is not in conflict with the data on the gravitational lensing. Then, 
the analysis of [18] supports the model considered in this paper, which includes general 
features of MOND, too. 

4. Conclusion 

Briefly summarizing the consideration, we have described the deflection of light by the dark 
mater halo modelled by the constant 4-vector fleld, in detail. The analytical results are in 
agreement with the limit of small velocity of rotation f q ^ in previous studies. 

The deflection angle for the inflnite halo is twice of the value for the global monopole. 
The effective lensing mass is factor 7r/2 greater than the dark matter mass at the same 
distance. For the typical rotation velocity in spiral galaxies, the characteristic value 
of deflection angle is about several arcseconds, that is consistent with values observed 
astronomically. The relative surface brightness of two images of lensed object is determined 
by the inverse ratio of impact angles with respect to the dark lens. 

The effect of halo edge is reduced to slow dependence of deflection angle on the impact 
distance, that weakly perturbs the gravitational lensing. 

In conclusion, we have to make several important comments: 

First, it would be wrong to identify the model of halo accepted in the paper with known 
model of isothermal halo. Such the opinion could not be correct because of the following 
reasons: 

i) The physics of vector fleld model in the paper and the physics of general isothermal 
model are certainly different. 

ii) The vector fleld model of dark matter includes inner and outer edges, so that between 
the edges the metric is described by the same metric as for the isothermal halo. Hence, 
the deflection at impact distances greater than the radius of inner edge and much less 
that the radius of outer edge coincides with the result of isothermal halo. That is 
a new statement for the limit of infinite halo in the vector field model with different 
physics. The final equation for the defiection is known [19], but the physical motivation 
is different. For instance, the model does not contain the singularity of isothermal halo 
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(that is quite evident). So, the identification with the singular isothermal halo looks 
misleading. 

iii) The effects related with the inner edge absent in the case of isothermal halo were 
investigated in the literature [18]. 

iv) The effects related with the outer edge absent in the case of isothermal halo as was 
represented in the literature [12], contain some errors, i.e. the presence of fracture 
point in the bending curve at the edge was missed. The representation of correction 
in the paper has required further application of method, which gives the correct result 
for the infinite halo, that has been usefully made in the initial part of Section 2. 

Thus, the paper describes the physics different from the singular isothermal halo. The 
results concerning the limit of infinite halo in the vector field model and the effects of outer 
edge of halo are new. 

Second, we have to emphasize that the reason for the difference between the effective 
mass deflecting the light at the impact distance r and the dark matter mass inside the 
same radius r by a factor found in the paper is quite clear: the deflection takes place due to 
the mass inside the radius r and effectively due to the dark matter mass outside the same 
radius. This reason is evident and transparent. 

This work is partially supported by the grant of the president of Russian Federation for 
scientific schools NSc-1303.2003.2, and the Russian Foundation for Basic Research, grant 
04-02-17530. 
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